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On a Si(lll) vicinal face near the structural transition temperature, the lxl structure and 
the 7x7 structure coexist in a terrace: the lxl structure is in the lower side of the step edge 
and the 7x7 structure in the upper side. The diffusion coefficient of adatoms is different in the 
two structures. Taking account of the gap in the diffusion coefficient at the step, we study the 
possibility of step wandering induced by drift of adatoms. A linear stability analysis shows that 
the step wandering always occurs with step-down drift if the diffusion coefficient has a gap at the 
step. Formation of straight grooves by the step wandering is expected from a nonlinear analysis. 
The stability analysis also shows that step bunching occurs irrespective of the drift direction if the 
diffusion in the lower side of the step is faster. The step bunching disturbs the formation of grooves. 
If step-step repulsion is strong, however, the step bunching is suppressed and the straight grooves 
appear. Monte Carlo simulation confirms these predictions. 

PACS numbers: 81.10.Aj, 05.70.Ln, 47.20.Hw, 68.35.Fx 
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I. INTRODUCTION 

On vicinal faces of Si(lll) 00,0 and Si(001) @, step 
wandering occurs at high temperatures when a specimen 
is heated by direct electric current. The current direction 
to cause the step wandering is step-down on the Si(lll) 
vicinal face 0, 0, and step-up on the Si(001) vicinal 
face 0. 

The cause of the step wandering is drift of adatoms by 
the current 0, |E 0J3 • The drift is in the same direction 
as the current 0, fjj- If the step is impermeable 0, , 
the step wandering occurs with step-down drift, as in the 
Si(lll) vicinal face [1 H 0. If there is alternation of 
the anisotropy in the diffusion coefficient on consecutive 
terraces flj|, as in the Si(OOl) vicinal face Q the step 
wandering occurs with step-up drift. 

On Si(lll) surfaces, the lxl structure is reconstructed 
and the 7x7 structure appears at low temperatures 
(< 860°C). In a vicinal face near the transition temper- 
ature, the 7x7 structure spreads from the upper side of 
the steps, and the two structures coexist in a terrace. Re- 
cently, Hibino and co-workers ffj observed step wander- 
ing near 860° C during growth. Due to the in-phase step 
wandering, grooves perpendicular to the steps appear on 
the vicinal face. Kato and co-workers 01 studied the 
step wandering theoretically. Focusing on the difference 
in diffusion coefficient of the two structures, they showed 
that the step wandering occurs in growth if the diffusion 
coefficient on the lxl structure is larger than that on 
the 7x7 structure. 

With the two phases coexisting, the drift of adatoms 
may also cause the step wandering instability on the 
Si(lll) vicinal face. In this paper, we study the pos- 



sibility of morphological instabilities induced by the drift 
of adatoms with the gap in the diffusion coefficient on 
the upper and the lower terraces. 



II. MODEL 

We consider a vicinal face where steps are running par- 
allel to the x-direction bordering terraces of a width I on 
average. The y-direction is chosen toward the step-down 
direction. If impingement and evaporation of adatoms 
are neglected, the adatom density c(r, t) is determined 
by 



dc{r,t) 
dt 



V-i(r,t), 



(1) 



where j(r, t) is the adatom current on the surface. With 
step-down drift, the adatom current is given by 

j(r,t) = -D s (r) (Vc(r,i) - ^ , (2 ) 

where D s (r) is the local diffusion coefficient, F the force 
to cause the drift and is positive for the step-down drift, 
and e. y the unit vector toward the step-down direction. 
We assume that the diffusion coefficient D s takes two 
values in a terrace: D s = D\ in the lower side of a step 

(n) 

edge, y n < y < y n + 1[ , and D s = D2 in the upper side, 
1 + < y < y n , where y n (x,t) is the position 



of the nth step and l^(x, t) is the terrace width of the 
lower side structure (1 x 1 in Si(lll)). 

Solidification and melting occur at step edges. In local 
equilibrium at a step, the adatom density is given by 
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(3) 



2 



where Cg q is the equilibrium adatom density at the iso- 
lated step, fl the atomic area, /3 the step stiffness, k the 
step curvature and U n the step-step interaction potential. 

From the continuity of the adatom current and the 
adatom density, the boundary conditions at the phase 
boundary are given by 



(4) 
(5) 



where fib is the normal vector of the boundary and +(— ) 
indicates the lower (upper) side of the step. 

By solving the diffusion equation Q in a static approx- 
imation with the boundary conditions, Eqs. the 
adatom density is determined. The normal step velocity 
V n is given by 

v n = m s -u\ Vn _-j\ yn+ ), (6) 

where n s is the normal vector of the step. 

III. STABILITY ANALYSIS 

When the steps and the boundaries of two phases 
are straight, the adatom density in the quasi-static ap- 
proximation is given by c(y) = Aq + Boe^^ v ~ Vn ' , where 
/ = F/k-gT. In the lower side of a step, the coefficients 
A and Bq are given by 



A 



B n = 



D 2 cUe^-l) 



D 2 e/* 2 ( e M n) - l) + D 1 ( e /4" ) - 1)' 



(D 1 - D 2 y> (e 



D 2 ef l 2 ) {e< } - I) + Di(eM" 5 - 1) : 



(7) 



(8) 



where = y n +i — Un and l 2 = — 4 ■ In the upper 
side of the step, the coefficients are given by the same 
form with the replacement D\ «-» D 2 and 4 <-> l 2 ■ 
The adatom current jo is constant on the whole of the 
n-th terrace and is given by 



Jo (I 



DiD 2 c°f(e 



fl [ 



1 



D 2 e^\e^ 



(9) 



l)+£>i(e^ - 1) 

Since we have neglected the step repulsion, the current 
is a function of and l 2 an d does not depend on the 
neighboring terrace widths. 

If the ratio of the width of the two structures 7 = 

4™' /4 i s nxC( i, Eq.© gives the current jo(l) as a func- 
tion of the terrace width. When the steps are equidis- 
tant, the adatom current at the step positions from 
the upper terrace equals to that onto the lower terrace, 
jo\ y = 3o\y n + - The velocity of the steps vanishes and 
the steps do not move. If the step interaction is neglected, 
the stability of the equidistant steps for step pairing is 
determined by since it controls the balance of the 



incoming and outgoing current with a pairing fluctua- 
tion of the terrace width (the repulsive interaction tends 
to stabilize the system). It is unstable for step pairing 
if j'o(0 > and stable otherwise. From Eq. ©, j' Q (l) is 
given by 



foQ) = 



(Di — D 2 )DiD 2 c° f 2 e fl2 hl 



eJ l - 1 



-1) 



■£>i(e^ 
I s 



I)] 2 



h 



(10) 



When the surface diffusion in the lower side of the step 
is faster than that in the upper side (Di > D 2 ), the 
vicinal face is unstable for the step pairing. The stability 
is independent of the ratio 7 of the widths and the drift 
direction. 
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FIG. 1: A terrace bounded by tilted steps with an angle 9. 

To find the condition for the in-phase step wandering, 
we consider an equidistant train of straight steps tilted 
with an angle 9 (Fig. 2J. The adatom current on the 

terrace in the x-direction, J x ^ is calculated as 

4°' = dy(j\\ cos6 ~ 3± sm8) 

Jy n 

(£»i - D 2 e f±h± )(e f±h± - l){e f±l2± - 1) 
£) 2e /'2( e /xiii - 1) + Di(ef^ - 1) 

x(£>i-£>2)c° q tan0, (11) 

where jy and j± indicate the adatom current in the £ and 
£ directions, f± = f cos9, l± = lcos9, li± = l\ cosQ and 
h± = l 2 cos8. When the step distance is small enough, 
f±l± <C 1, it becomes 



j(o) 



(D 1 - D 2 ) 2 fhl 2 c° eq 



{D x h + D 2 l 2 ) 1 + (r] x y 



(12) 



where we assume y n (x, t) = nl + rj{x, t) and rj x = di]/ dx. 
In addition to Eq. i|12|) . there is current Jx caused by 
the change of the chemical potential along the step: 



J« = - cos 2 6(D 1 l 1± + D 2 l 2± )c° c(1 ^- (-^ 



, (13) 
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where fi = Q[3k. 

The evolution of the step position is determined by the 
adatom current in the x-direction as 



dr) 
dt 



dU o) +jn 



ID- 



d_ 

dx 



dx 
a 2 rj x 



Q!4 3 



I + 772 l + rgdx V(l+ rgfl 2 
where the coefficients a 2 and 0:4 are 

,(£>!- D 2 ) 2 fhl2C° eq 



a 2 



0- 



(ZM 2 + D 2 h) 

^ 



a 4 = n{Dih + D 2 l 2 ) 



cq fc B T' 



If the step position is of the form rj{x) 
the linear amplification rate is 



a 2 q 



4 



,(14) 

(15) 
(16) 

(17) 



The coefficient 014 is always positive and suppresses the 
step fluctuation. With step-up drift (/ < 0) the coeffi- 
cient a 2 is negative and suppresses the step fluctuation, 
while a 2 becomes positive and the step wandering occurs 
with step-down drift (/ > 0). 

Equation (|14|) is the same type of equation descri bing 
the step wandering in other conserved systems [T3L Il4l 
IT5I Il6j . The solution of the equation shows a regular 
periodic pattern whose amplitude increases in a power 
law of time as t 1 / 2 [l4j|. As a result periodic grooves will 
be formed. 



IV. MONTE CARLO SIMULATION 

We perform Monte Carlo simulation for solid-on-solid 
steps of a square lattice model. The boundary condi- 
tion is helical in the y-direction and periodic in the x- 
direction. We assume that 7 is fixed to 1 so that the 
phase boundary is at (y n + y rl +i)/2 when the steps move. 
In the lower side of a step (supposedly the lxl region), 
an adatom on the site moves to (i ± with the 
probability pd = D i/4 and to (i , j ± 1 ) with the probabil- 
ity Pd = Di(l ± /a/2)/4. In the upper side of a step (the 
7x7 region), the parameter D\ is replaced by D 2 . The 
diffusion across the boundary of the two regions takes D 2 . 
In our simulations, we assume that (D\,D 2 ) = (l,a) or 
(a, 1) with a < 1. 

Solidification and melting occur at the lower edge of 
the step positions ■ The probabilities for solidification 
p + and melting p_ arc given by 



P± 



1 + cxp 



( AE S + AUt< 



(18) 



where A£" s is the increment of the step energy, <j) the po- 
tential gain by solidification and At/ is the change of the 



step-step interaction potential. We assume the repulsive 
interaction potential U n of the nth step takes the form 



A 



(19) 



where y n (xi) is the position of the nth step at x = a;,-. 

We first carry out the simulation with the diffusion 
coefficients (Di,D 2 ) = (1,0.1), the stiffness ft/k^T = 
1.64 and the equilibrium density c° q = 0.18. The system 
size is 512 x 512 and the initial step distance is I = 16 
(the step number is N = 32). Initially, the steps are 
straight and equidistant, and there are a few adatoms on 
the terraces. 
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FIG. 2: Snapshots of a destabilized vicinal face for D\ > D2: 
with (a) step-up drift and (b) step-down drift . 

When we neglect the step-step repulsive interac- 
tion, the vicinal face is unstable for the step bunching 
(Fig. HJa)), which agrees with the analysis in Sec. II I II 
The step bunching occurs irrespective of the drift direc- 
tion, but the form of the bunches changes with the drift 
direction. The bunches are straight with step-up drift 
(Fig-GI a )) an d wander with step-down drift (Fig. Gib)). 
With a strong repulsive interaction, A = 300, the step 
bunching is suppressed and the in-phase step wander- 
ing occurs. We have not studied the growth laws of the 
step width and the period of grooves, but the form of 
the groo ves (Fig. |3J) qualitatively agrees with the solu- 
tion [11 [13 of Eq. JTJ. 
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FIG. 3: Snapshot of step wandering induced by step-down 
drift with a strong step-step repulsion for D\ > D2- Except 
for the repulsive potential, the parameters are the same as in 

Fig. □>) . 
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FIG. 4: Snapshots of a vicinal face for D\ < D2 : with (a) step- 
up drift and (b) step-down drift. The repulsive interaction is 
absent. 



(Fig.|3Ja)), neither step wandering nor step bunching oc- 
curs. With step-down drift, the step wandering occurs, 
but no indication of bunching is seen (Fig. Efb)). As 
time passes, periodicity selection proceeds, resulting in 
straight grooves as those of Fig. 



V. SUMMARY AND DISCUSSION 

In this paper, we studied the drift-induced morpholog- 
ical instabilities on a vicinal face with two phases. As 
is seen from Eq. i|15|> . the step wandering occurs with 
step-down drift unless D\ — D^. If Di = D2, our model 
reduces to a model of a simple vicinal face. We have 
already found |lOj | that, with step-down drift, wander- 
ing instability occurs in such a simple vicinal face if the 
kinetic coefficient of the step is finite. Although in the 
present paper we have assumed an infinite kinetic coef- 
ficient (local equilibrium at the steps), it must be finite 
in reality. Therefore we may say wandering instability is 
expected irrespective of the diffusion ratio. In all cases 
straight grooves parallel to the drift is produced due to 
the in-phase step wandering in Monte Carlo simulation 
if step bunching is suppressed. 

In the Si(lll) vicinal face near the transition temper- 
ature, the diffusion coefficient in the lower side of a step 
is larger than that in the upper side [l^. With the as- 
sumption that the ratio of the widths of the two struc- 
tures is constant, the step bunching occurs irrespective 
of the current direction and the step wandering occurs 
with step-down current. In a vicinal face of large incli- 
nation, however, the step bunching is suppressed due to 
the strong step-step repulsion, and only grooves induced 
by the step wandering may be observed with step-down 
current. 

In our model, we assumed that the boundary of the 
two structures moves in concert with the steps and the 
ratio of the widths of two structures is constant. The 
drift direction to cause the instabilities does not depend 
on the ratio of the widths. In reality the ratio of the 
widths changes with temperature and the motion of 
the boundary does not automatically follow the steps. 
To study the morphological development of this system 
in detail, we need to extend our model to include the 
freedom of the motion of the boundary [2(J • 



Acknowledgments 



We also carry out the simulation with D\ < D2 
(Fig.0}. The diffusion coefficients are {D ll D 2 ) = (0.1, 1) 
and the step-step repulsion is neglected (A = 0). Other 
parameters are the same as Fig. With step-up drift 
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